Astronomy & Astrophysics manuscript no. DelVel'astroph 


©ESO 2011 


June 14,2011 





in 



Oh- 



oo 



o 



X 



The nearby eclipsing stellar system 6 Velorum 
III. Self-consistent fundamental parameters and distance* 

A. Merand 1 , P. Kervella 2 , T. Pribulla 3 - 4 , M. G. Petr-Gotzens 5 , M. Benisty 6 , A. Natta 7 ' 8 , G. Duvert 9 , D. Schertl 10 , and 

M. Vannier 11 



European Southern Observatory, Alonso de Cordova 3107, Casilla 19001, Vitacura, Santiago 19, Chile 

2 LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Universite Paris Diderot, 5 place Jules Janssen, 92195 Meudon, France 

3 Astronomical Institute Slovak Academy of Sciences 059 60 Tatranska Lomnica, Slovak Republic 

4 Astrophysikalisches Institut und Universitats-Sternwarte, SchillergaGchen 2-3, 07745 Jena, Germany 

5 European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany 

6 Max-Planck-Institut fur Astronomie, Konigstuhl 17, D-69117 Heidelberg, Germany 
s ! , 7 Osservatorio Astrofisico di Arcetri, INAF, Largo E.Fermi 5, 1-50125 Firenze, Italy 

f— | ■ 8 School of Cosmic Physics, Dublin Institute for Advanced Studies, Dublin 2, Republic of Ireland 

9 UJF-Grenoble 1/CNRS-INSU, Institut de Planetologie et d'Astrophysique de Grenoble (IPAG) UMR 5274, Grenoble, F-38041, 

^— . ' France 



Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

Laboratoire Fizeau, Universite de Nice, CNRS-Observatoire de la Cote d' Azur, 06108 Nice Cedex 2, France 



Received — ; accepted — 

ABSTRACT 



Context. The triple stellar system S Vel (composed of two A-type and one F-type main sequence stars) is particularly interesting as it 
contains one of the nearest and brightest eclipsing binaries. It therefore presents a unique opportunity to determine independently the 
physical properties of the three components of the system, as well as its distance. 
Aims. We aim at determining the fundamental parameters (masses, radii, luminosities, rotational velocities) of the three components 
5— l ■ of S Vel, as well as the parallax of the system, independently from the existing Hipparcos measurement. 

^ ' Methods. We determined dynamical masses from high-precision astrometry of the orbits of Aab-B and Aa-Ab using adaptive optics 

(VLT/NACO) and optical interferometry (VLTI/ AMBER). The main component is an eclipsing binary composed of two early A- 
type stars in rapid rotation. We modeled the photometric and radial velocity measurements of the eclipsing pair Aa-Ab using a self 
consistent method based on physical parameters (mass, radius, luminosity, rotational velocity). 

Results. From our self-consistent modeling of the primary and secondary components of the 6 Vel A eclipsing pair, we derive their 
fundamental parameters with a typical accuracy of 1%. We find that they have similar masses, respectively 2.43 ± 0.02 M Q and 
2.27 + 0.02 M Q . The physical parameters of the tertiary component (S Vel B) are also estimated, although to a lower accuracy. We 
obtain a parallax n = 39.8±0.4mas for the system, in satisfactory agreement (-1.2cr) with the Hipparcos value (^Hi P = 40.5 ±0.4 mas). 
Conclusions. The physical parameters we derive represent a consistent set of constraints for the evolutionary modeling of this system. 
The agreement of the parallax we measure with the Hipparcos value to a 1% accuracy is also an interesting confirmation of the true 
. accuracy of these two independent measurements. 



Key words. Stars: individual: (HD 74956, 6 Vel); Stars: binaries: eclipsing; Stars: early-type; Stars: Rotation; Techniques: high 
angular resolution; Techniques: interferometric 

1 . Introduction pie is provided by the recent interferometric observations of the 

A0V benchmark star Vega, that confirmed that Vega, as previ- 



^ . Early-type main sequence stars exhibit a number of peculiar- ously shown by Gullivel - Hill & Adelman ([19941 is a pole 

ities usually not encountered in cooler stars: fast rotation, de- on fast rotator near critical velocity (Aufdenberg et aL 2 006). 

bns disks, enhanced surface metalhcities (Am), magnetic fields The same i nter ferometric observations showed that Vega har- 

and rapid oscillations (Ap and roAp stars), etc. Although stellar borfj a hot debris disk whhin within 8 AU from the star (Abs0 

structure and evolution models are now rather successful in re- et a j | 2QQ6l i 
producing the observed physical properties of most A-type stars, 

the observational constraints on these models remain relatively 8 Vel ( HD 74956, HIP41913, GJ 321.3, GJ 9278) is abright 

weak, occasionally leading to surprising discoveries. An exam- multiple star including at least three identified components, and 

is among our closest stellar neighbors, with a revised Hipparcos 

Send offprint requests to: A. Merand parallax of 7r Hi p = 40.5 + 0.4 mas (van Leeuwen |2007]l. This 

* Based on observations made with ESO telescopes at Paranal ob J ect has man y observational peculiarities. Firstly, it was dis- 

Observatory, under ESO programs 076.D-0782(B), 081.D-0109(B), covered only in 1997 that 6 Vel hosts one of the brightest of all 

081.D-0109(C), 282.D-5006(A) and Arcetri GTO program 084.C- known eclipsing binaries (Otero et al. l20001 >. with a remarkably 

0170(C) long orbital period (P « 45 days). This eclipsing binary is also 

Correspondence to: amerand@eso.org one of the very few that are easily observable with the naked eye 
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(my « 2). The eclipsing pair was first resolved using optical in- 
terferometry by Kellerer et al. (I2007I I. Secondly, 5 Vel is known 
to have a moderate thermal infrared excess (e.g. Aumann 1985 
Su et al. 2006), and Spitzer observations revealed a spectacu- 
lar bow shock caused by the motion of 6 Vel in a dense inter- 
stellar cloud (Gaspar et al. |2008 ). The presence of interstellar 
material was also reported by Hempel & Schmitt (2003 ), who 
observed two red-shifted absorbing components in absorption in 
the Ca II K line, of probable interstellar origin. In Paper I of 
the present series, Kervella et al. (2009) confirmed that the in- 
frared excess is essentially emitted by the bow shock, and not 
warm circumstellar material located close to the stars. In the 
framework of a search for resolved emission due to debris disks, 
Moerchen et al. (120101 ) obtained thermal infrared images of S Vel 
using the Gemini South telescope and the T-ReCS instrument, 
and detected a marginally resolved emission at A — 10.4 pm. 

In Paper II, Pribulla et al. (201 1 ) used a combination of high- 
resolution spectroscopy and photometric observations (from the 
SMEI instrument, attached to the Coriolis satellite) to derive an 
accurate orbital solution for the eclipsing binary 6 Vel A, and 
estimate the physical parameters of 5 Vel Aa and Ab. They iden- 
tified that the two eclipsing components are fast rotating stars, 
with respective masses of 2.53 + 0.11 M Q and 2.37 + 0.10 M 
(* 4% accuracy), and estimated the mass of <5Vel B to be 
~1.5M a . 

In spite of this recent progress, uncertainties remain on the 
fundamental parameters of the different components of the sys- 
tem, in particular on their exact masses. Taking advantage of 
the availability of NACO astrometry of the 5 Vel A-B pair, and 
new interferometric observations from the VLTI/AMBER in- 
strument, we propose here to revisit the system along two direc- 
tions. In Sect. [2] we describe our new VLTI/AMBER interfero- 
metric data, as well as our re-analysis of the spectroscopic and 
photometric data previously used in Paper II. Sect. 12.21 is ded- 
icated to the description of our self-consistent model, and the 
derivation of an improved orbital solution, physical parameters 
as well as an independent distance. In Sect. |3]we employ NACO 
astrometry of the visual 5 Vel A-B binary to obtain an improved 
orbital solution. Compared to our work presented in Paper II, 
this new analysis result is a clearer view and better confidence in 
the derived fundamental parameters of the system (for all three 
components), thanks to the redundant nature of our data and our 
independent determination of the distance. 

2. The orbit and parameters of 5 Vel Aa and Ab 

2.1. Observations and data analysis 
2.1.1. Interferometry 

AMBER (Petrov et al. I2007I ), the three-telescope beam com- 
biner of the VLTI, has the proper angular resolution to resolve 
the Aa-Ab pair. This instrument combines simultaneously 3 ATs 
(Auxiliary Telescopes) or 3 UTs (Unit Telescopes) of the VLT 
and operates in the near infrared (H and K band). It has a choice 
of spectral resolutions of R ~ 35, R ~ 1500 or R ~ 15000. 
For this study, we had data in low resolution (H+K bands at 
R ~ 35) and medium resolution (H or K band, at R ~ 1500). We 
used baselines of the order of 100 m in order to obtain spatial 
resolution in the milli-arcsecond regime. These interferometric 
data have been collected in a dedicated program (ESO program 
076.D-0782), as well as during Guaranteed Time (GTO) from 
Arcetri Observatory. We present here a reduction of these data. 

We reduced the data using the AMBER reduction package 
amdlib3 (Chelli et al. 120091 Tatulli et al. 120071 and performed 



Table 1. AMBER separation vectors (primary to secondary) 



date 


toward East 


toward North 


MJD 


mas 


mas 


53427.09 


-4.1+0.4 


7.3 ±0.5 


53784.08 


-6.5+0.5 


13.8 ±0.3 


54819.25 


-7.3 ±0.2 


16.1 ±0.2 


54832.12 


2.7 + 0.7 


-6.4 ± 0.4 


55147.38 


1.6 + 0.5 


-4.0 ± 0.5 



Notes. For simplicity, error bars have been estimated along E and N 
directions even though the true error is an ellipse with orientation dieted 
by the x 1 map. It turns out the ellipses have small flattening, so the 
approximation is relevant. 



the calibration using stellar calibrators chosen in the catalogue 
by Merand et al. (2005) and a custom software which estimates 
and interpolates the transfer function of the instrument. For each 
night, we derived the separation of Aa and Ab using a.y 2 map as 
a function of the separation vector (two parameters). The other 
parameters, such as flux ratio or individual diameters, were set 
using simple hypothesis and their choice did not affect signifi- 
cantly our final estimated angular separations. The resulting sep- 
aration vectors are listed in Table Q] in coordinate towards East 
and North (which correspond to the u and v axes in the projected 
baselines map). The error bars on this vectors were estimated in 
the^ 2 map. 

2.1.2. Photometry 

We used the photometric data from the SMEI satellite, presented 
in Paper II. The available quantity is the relative flux normalized 
to the value outside of the eclipses, since there was no absolute 
calibration of SMEI data available. We corrected for the pres- 
ence of the B component which is in the field of view of SMEI. 
From our model of B, the expected flux ratio between B and Aab 
is 7.5% in the SMEI bandpass. The transmission of the instru- 
ment has a triangular shape that peaks at 700 nm, with a quantum 
efficiency around 47%, and falls to « 5% at 430 nm towards the 
blue, and 1025 nm towards the red (Spreckley & Stevens 2008). 
We removed this contribution which, if not taken into account, 
would result in an underestimation of the depth of the eclipses. 
We also incorporated to our photometric dataset the photomet- 
ric measurement we derived of Aab in the K band (Paper I). We 
use this value in our fit as the only constraint in term of absolute 
photometry. 

2.1.3. Spectroscopy 

The observables we derived from the visible spectroscopic data 
are the broadening functions (BF) presented in paper II (see this 
reference for more explanations). These functions contain a lot 
of information: not only do they contain the radial velocities that 
result from the orbital motion, but also the broadening due to the 
stellar rotation and the flux ratio in the considered band. 

From the observed BF, it is possible to derive the v sin i from 
the two components. After a few experimentations using the stel- 
lar surface model we are going to present later, we found that the 
following ad-hoc function parameterizes well the BF for a star 
seen from the equator: 



BF(v) 
BF(v ) 
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0.5 

orbital phase 

Fig. 1. 6 Vel Aa-Ab spectroscopic broadening functions (BF). 
Upper: BF as a function of orbital phase; Lower: residuals after 
the fit. See Fig. [2] for the fit at the phase corresponding to the 
dotted line (<p «s 0.17). Our method to extract the radial velocity 
does not take into account the eclipses, as can be seen by the 
large residuals during the eclipses (0 w 0.40 and <p & 0.97, upper 
right panel). 



where a is the only parameter constraining the gravity dark- 
ening and Vo is the velocity offset. The function is defined for 
|v — vol < v sin i only, and its value is otherwise. Using this ana- 
lytical model and a global fit, we estimated v sin i for each com- 
ponent and the radial velocities for each epoch (see Fig.[T]and|2] 
for the quality of the fit). We found v sin i to be 143.5 + 0.2 km/s 
and 149.6 + 0.2 km/s for Aa and Ab respectively. Incidentally, 
we find a a = 0.460 + 0.003 and a b = 0.45 1 + 0.003. The rotation 
rate value is relatively independent of the actual gravity darken- 
ing (parameterized here by a) since it is set by the width of the 
broadening function, not its shape. 

For our fit of the orbit and the stellar parameters, we do 
not use the center-to-limb darkening we derive here from the 
broadening functions. The surface brightness distribution is con- 
strained by the photometric profil of the eclipses. However, we 
will check a posteriori the agreement between our best fit model 
and the limb-darkening derived in the analytical BF by modeling 
the BF from our model. See Sect. 12.2.31 and, more specifically, 

Fig.m 

2.2. Global fit 

2.2.1 . Self consistent model 

In order to extract the fundamental parameters (masses, radii, 
surface temperatures, semi-major axis, etc.) from the obser- 
vational data, we propose a self consistent modeling centered 
around the use of physical quantities: we model the system us- 
ing two stars whose characteristics are computed based on their 
radii, total luminosities and mass. 




Vrad (km/s) 

Fig. 2. Example of one fit of a broadening function (for the 
phase represented as a dotted line in Fig. \T\ <t> « 0.17). The 
analytical fit uses the sum of two functions described in eq. Q] 



To illustrate the advantage of this approach, we can con- 
sider that in order to model the eclipses, we could use an ad-hoc 
model based on fractional radii (ratio to the semi-major axis) and 
brightness ratio, but this would not lead directly to the funda- 
mental parameters of the system such as effective temperatures 
or luminosities. Our approach uses radii, masses and luminosi- 
ties: we get the fractional radii by self consistency between the 
semi-major axis based on Newton's form of the Kepler's law 
(from the masses and the period of the orbit) and the measured 
apparent semi-major axis (constrained by interferometric sepa- 
ration vectors). The brightness ratio arises from the luminosity 
and radii, and the photospheric models we use to model the sur- 
face of the stars. 

Our stellar surface model also includes stellar rotation. We 
model the appearance of the star using a model developed 
contemporarily and similar to the one used in Aufdenberg et 
al. ((2006) to model the interferometric visibilities of the star 
Vega. To compute the photometry, in particular during the 
eclipses, we generated synthetic images and integrated them to 
derive the light curves. The parameters we use are: 

- total mass of the system (Aa+Ab); 

- the fractional mass of Aa to the total mass of the eclipsing 
system (Aa+Ab); 

- the physical radii of each component; 

- the absolute luminosity of each component; 

- the v sin i of each component, to parameterize the rotation; 

In addition, we have the usual 7 parameters for the visual orbit: 

- the period; 

- the date of passage at periastron; 

- the eccentricity; 

- three angles: inclination, a> and Q; 

- the apparent semi-major axis (in milliarcseconds) 

The only parameterization of the physics of the two com- 
ponents is contained in the luminosity: the model we use of 
the stellar surface is the Roche approximation, which only uses 
the mass, radius, luminosity and rotational velocity. Once the 
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Table 3. Derived fundamental parameters for S Vel Aa and Ab.. 
Parameters for Vega are displayed for comparison. 





unit 


6 Vel Aa 


6 Vel Ab 


Vega 


Mass 


M e 


2.43 ± 0.02 


2.27 ± 0.02 


2.3 ± 0.2 


Luminosity 


U 


67 + 3 


51+2 


37 + 3 


Polar Radius 


l< 


2.79 ± 0.04 


2.37 ± 0.02 


2.26 ± 0.07 


Equ. Radius 


Re 


2.97 ± 0.02 


2.52 ± 0.03 


2.78 ± 0.02 


Polar r cff 


K 


10100 


10120 


10150 


Equ. 7W. 


K 


9700 


9560 


7900 


Avg. r eff . 


K 


9450 


9830 


9100 






0.61 


0.60 


0.91 


Polar log(g) 


cm s~ 2 


3.90 


4.10 


4.10 + 0.1 


Eq. log(g) 


cm s~ 2 


3.78 


3.99 


3.65+0.1 


/ 


deg 


~ 90 


~ 90 


4.7 ± 0.3 


rotation rate 


1/d 


0.95 


1.17 


1.90 


metallicity 


[M/H] 


-0.33 1 


-0.33 1 


-0.5 1 



Notes. Parameters for Vega are adapted from Aufdenberg et al. ( 2006 ). 
The references for the metallicity are 'Gray et al. ( 1200 6 1 and 2 Castelli 
& Kurucz ( [19941 ). 



shape of the surface is computed, we link the local surface grav- 
ity and local effective temperature using the von Zeipel the- 
ory (von Zeipel 119241 and Aufdenberg et al. 120061 for more de- 
tails). The luminosity is constrained using several mechanisms: 
through the absolute photometry of the system, but also through 
the surface brightness that sets the depth of the eclipses. 

The advantages of using the apparent semi-major axis com- 
pared to the physical quantity are simple. First of all, it is di- 
rectly related to one of our observables: the interferometric sep- 
aration vector. Secondly, we already have the semi-major axis 
by combining the Kepler's third law, the period and the total 
mass. Combined with the angular semi-major axis, we derive a 
distance independently from the Hipparcos value (a brief discus- 
sion is presented in Sect. I2.31 >. The distance is derived internally 
in our model and used to extract a model apparent magnitude in 
the K band, which is used as one of the constraints as we men- 
tioned before. 

It is to be noted that we make the following assumptions: 

- the stars have their rotation axis perpendicular to the plane 
of the orbit. 

- we use a Von Zeipel gravity darkening coefficient of /3 = 0.25 
in T e g oc g^ ff , where T e g and g e g are the effective tempera- 
ture and gravity at the surface. Even if recent observational 
constraints from Monnier et al. (120071 1 suggest f3 ~ 0.19, we 
chose to use von-Zeipel's classical value since it does not 
lead to qualitatively nor quantitatively different results in our 
case. 

2.2.2. Orbital and stellar derived parameters 

The result of the global fit is in excellent agreement with all the 
observed data we used for the fit: the relative photometry of the 
eclipses (Fig.[3]and|4|, the radial velocities (Fig. [5]) and the sepa- 
ration vectors (Fig.[6]l. The corresponding orbital parameters and 
stellar parameters are presented in Tables [2] and [3] 

2.2.3. a posteriori verifications 

Our model is rather simple and does not take into account one 
aspect: the heating of one star by the other's radiation. This ef- 
fect could be important in our case in different conditions. The 
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0.98 



0.99 



Fig. 3. 6 Vel Aa-Ab SMEI photometric measurements (points) 
and fit of the eclipses (line) as a function of orbital phase: the 
upper panel is the the primary eclipse, the lower panel the sec- 
ondary. 



first one is if one star was heated by the other one and devel- 
oped a bright spot on its surface. We can discard this possibil- 
ity because of the absence of photometric variations outside the 
eclipses. The second case where the heating by the other star can 
be a problem is if this contribution is enough to actually mod- 
ify the temperature structure of the stellar photosphere. We can 
check a posteriori that the radiation received from the other star 
is of the order of 3% of the radiation emitted (assuming similar 
surface brightness and a 2.5 R e star seen at a distance of 89 R Q ). 
Approximating the surface as a blackbody, it corresponds to an 
increase of temperature of less than 1% (1.03^ 4 ). We can thus 
assume that our hypothesis that locally the photosphere can be 
approximated by a ATLAS model is correct. 

We can also check the consistency of our model beyond 
the fit of the data we presented. In particular, we can compare 
the predicted broadening functions based on our model and the 
broadening functions we observe since we did not implement the 
direct fit of broadening function to our model. Doing so (Pig. [7j, 
the comparison is very satisfactory, even though the wings of the 
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Table 2. Aa-Ab orbital and stellar parameters fitted on the photometric, interferometric and radial velocity data, using our model. 





Aa Ab 


constrained from 


a (mas) 


16.51 ±0.16 


interferometry 


Total mass (M Q ) 


4.69 ± 0.03 


all observations 


M A a/(M A a + M Ab ) 


0.516 ±0.001 


spectroscopy 


Polar Radius (R ) 


2.79 ± 0.04 2.37 ± 0.02 


photometry 


Luminosity (L Q ) 


67 ±1 51 ±1 


photometry 


v sin i (km/s) 


143.5 ±0.2 149.6 ±0.2 


spectroscopy 


Period P (d) 


45.1503 ±0.0002 


all observations 


MJD modulo P 


19.159 ± 0.010 


all observations 


e 


0.287 ± 0.001 


all observations 


i (deg) 


89.04 ± 0.03 


all observations 


0) (deg) 


109.79 ± 0.09 


all observations 


a (deg) 


65.0 ± 0.6 


interferometry 


V y (km/s) 


-9.78 ± 0.07 


spectroscopy 


X 1 Kad Primary (err 0.75km/s) 


1.13 




X 1 Kad Secondary (err 1.3km/s) 


1.03 




X 1 Photometry (err 0.75%) 


1.15 




X 2 Interferometry 


1.05 





Notes. The last four lines present the agreement as reduced^ 2 . 




Fig. 4. 6 Vel Aa-Ab synthetic images in linear scale, close to the eclipses: the upper panel is the configuration of the primary eclipse, 
the lower panel is the secondary eclipse. 
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orbital phase 

Fig. 5. 5 Vel Aa-Ab radial velocities as a function of orbital 
phase, with model from the fit over plotted and residuals plot- 
ted in the lower panel. Round points (blue) are for the primary 
and diamonds (red) for the secondary. 




mas 

Fig. 6. 5 Vel Aa-Ab visual orbit: AMBER positions (black dots 
with error bars, see Table [TJ; fitted orbit (gray line), Aa stellar 
disk (dark gray disk with blue surrounding), the position of Ab 
at periastron (red dot) and Ab stellar disks for the positions of 
the eclipses (gray disk with red surroundings). Residuals to the 
orbit (red lines, too small to see in most cases) and quality of the 
fit (in number of sigmas). The overall agreement corresponds to 
a reduced^ 2 of 1.05. 



data (represented by the analytical function fit to the data in thick 
gray line) seems deeper than for the model. In other words, the 
gravity darkening of the model is slightly underestimated, but 
only by a small fraction, considering that a model without grav- 
ity darkening (small gray dots on Fig. O produces a very strong 
disagreement. 

Investigating the possible causes of this problem, we realized 
that we can reproduce the more pronounced darkening of the 
equator compared to the pole by tilting the star to 10 degrees 
from the plane of the sky: this makes the pole more in line of 
sight of the observer and hence increases the contrast between 
the pole and the equator. 

Forcing the inclination of the spin of the stars to be 80 de- 
grees instead of 90 degrees does not change dramatically the 
fundamental parameters of the star estimated from our fit. One 
of the reasons is that it changes the sin(i) by only 1.5%, the ac- 
tual rotational velocity is mostly unaffected. The fit converges as 
well as in the case of aligned spins, with fundamental parameters 
within error bars of the one estimated in the case we presented 
in the main part of this work. 

In conclusion of the analysis of the broadening function pre- 
diction of our model compared to the observed one, we find 
a confirmation of the consistency of our model with the data. 
This comparison may also indicate that our model underestimate 
slightly the gravity darkening, or, alternatively, that the stars 
have their rotational axis tilted on the order of 10 degrees toward 
the observer, which does not impact qualitatively nor quantita- 
tively the fit of the data that led to the estimation of the funda- 
mental parameters presented earlier. 

Another test is to compare the predicted wide band photom- 
etry with the observed ones. In the literature, there is a handful 
of data available, for the combined AB or A and B separately. 
Doing so, we see (Table |4|i that the bluest magnitudes are not 
reproduced well. Assuming aB -V excess of E(B - V) — 0.055, 
the difference is nicely explained using an ISM extinction law as 
presented in Kervella et al. (120041 1 . 

2.3. Distance 

Our parameterization allows us to derive the distance as the ra- 
tio between of the semi-major axis of the eclipsing component 
(from the total mass, the period and Kepler's third law) and 
its apparent semi-major axis (from the interferometric observa- 
tions). Such a distance estimate is particularly interesting as it is 
purely geometrical, and independent of the Hipparcos measure- 
ment. From our model, we obtain a parallax of n = 39.8+0.4 mas 
for 5 Vel. 

Comparing this value to the ^eip = 40.5 + 0.4 mas revised 
Hipparcos parallaxQ obtained by van Leeuwen (2007) shows a 
good agreement of the two measurements, within 1.2 cr. This 
confirmation of the true accuracy of these independent mea- 
surements, at the 1 % level, shows that the Hipparcos measure- 
ment was not disturbed by the binary nature of 5 Vel A. This 
somewhat surprising result is due to the similar brightness ratio 
LAa/^Ab ~ 1 .3 and mass ratio Mp^/Mp^ ~ 1.1 of the 5 Vel A 
pair. This results in a very small apparent displacement of the 
center of light of the Aab system during the orbit, with respect 
to the center of gravity of the two stars. Using our model of the 
eclipsing system, we computed the expected photocenter dis- 



1 The revised Hipparcos parallax is consistent within 1 cr with the 
original Hipparcos reduction (40.9 ± 0.4 mas; ESA 1997), and with 
the ground-based parallax of this star (49.8 ± 9.4 mas; Van Altena et 
al. fT995l 
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Table 4. Computed and observed broad band photometry. 









model (observation) 




A 


extinction 


band 


Aa 


Ab 


Aab 


B 


AB 


obs.-mod. 


E(B -V) = 0.' 


B 


2.39 


2.71 


1.78 


6.10 


1.76 (2.00 1 ) 


0.24 


0.23 


V 


2.41 


2.73 


1.81 


5.59 (5.54 2 ) 


1.78 (1.96 1 ) 


0.18 


0.18 


J 


2.43 


2.76 


1.83 


4.72 


1.76 (i.77 3 ) 


0.02 


0.04 


H 


2.45 


2.78 


1.85 


4.44 


1.75 (1.76 3 ) 


0.01 


0.03 


K 


2.46 


2.78 


1.85 (1.86 1 ) 


4.42 (4.40 2 ) 


1.76 (1.72 3 ) 


-0.04 


0.02 



Notes. The magnitudes predicted by the best fit model are listed together with the observed magnitudes (brackets). In bold is the K magnitude of 
Aab which was used to constrain our fit. The B column comes from the fit presented in Sect. [3] The last two columns are the observed differences 
and the expected extinction for E(B - V) = 0.055, which best fits the observed colors. Observed magnitudes are from Morel & Magnenat (1978 1 
1 , Paper I 2 and Skrutskie et al. ( 12006b 3 



placement during an full orbit. We find that the peak-to-peak: 
photocenter displacement is of the order of one milliarcsecond, 
which is much smaller than the apparent astrometric shift due 
to the parallax. The binarity of the system therefore did not bias 
significantly the Hipparcos parallax measurement, neither did 
the low brightness of the B component. The observations of the 
photocenter displacement through high-precision differential as- 
trometry with the VLT/NACO instrument will be the subject of 
a future article. 



3. The orbit and parameters of 8 Vel B 

3.1. Astrometric data 

The binarity of <5 Vel was discovered by S. I. Bailey in 1894 
from Arequipa, Peru (and independently by Innes 1895). Over 
more than one century, the separation between 6 Vel A and B has 
been decreasing at a rate which nicely matches the progression 
of the angular resolution of the successive generations of imag- 
ing instruments (visual observations, photography, electronic de- 
vices). This progression allowed a relatively regular tracing of 
the visual orbit of the pair, down to the sub-arcsecond separa- 
tions that occur around the periastron passage. With the advent 
of speckle interferometry (Tango et al. 1979 ) and the Hipparcos 
satellite (ESA 1997) the accuracy of the measured relative po- 
sitions improved significantly. In Paper I, we present in details 
the new data we obtained with the Very Large Telescope, using 
both the K-band adaptive optic system VLT/NACO (Rousset et 
al. l2003l Lenz en eta l. [T99"8j l and the N-band camera VLT/VISIR 
(Lagage et al. 2004). Thanks to the large aperture of the tele- 
scopes and the diffraction-limited angular resolution, these ob- 
servations provide us with new high-precision astrometry of the 
A-B pair. The resulting separations of S Vel B relatively to A are 
presented in Table [5] For the conversion of the separation mea- 
sured in pixels to angular separations, we adopted the pixel scale 
of 13.26 ± 0.03 mas/pixel (Masciadri et al. EOTjBl for NACO and 
75 + 1 mas/pixel for VISIR. The assumed NACO plate scale is in 
good agreement with the calibration by Neuhauser et al. (2008 ), 
who demonstrated that this figure is stable over a period of at 
least 3 years. The VISIR plate scale uncertainty is set arbitrarily 
to s 1%, although it is probably better in reality. The angular 
separation was only « 0.6" for the epoch of our observations. 

In addition to these new astrometric measurements, we also 
take advantage of the historical astrometric positions assembled 
by Argyle et al. (2002 ) in his Table 5, that includes 17 epochs be- 
tween 1895 and 1999. It is to be noted that these authors used the 
two speckle interferometry epochs from Tango et al. ( 1979 ) with 



Table 5. Differential astrometry of 5 Vel from NACO and VISIR 
images. 



UT date 


MJD-54 000 


[a B - a A ] 


[8 b- 


S A ] 






mas 




mas 


2008-04-01" 


557.0224 


-430.6 ± 


1.0 


457.8 


± 1.0 


2008-04-04" 


560.9976 


-430.7 ± 


1.0 


457.1 


± 1.0 


2008-04-06" 


562.0121 


-430.8 ± 


1.0 


457.1 


± 1.0 


2008-04-07" 


563.0048 


-430.7 ± 


1.0 


457.0 


± 1.0 


2008-04-20" 


576.9715 


-430.3 ± 


1.0 


452.9 


± 1.0 


2008-04-23" 


579.0231 


-430.0 ± 


1.0 


452.0 


± 1.0 


2008-04-24 1 ' 


580.0503 


-425.9 ± 


4.3 


454.0 


±4.5 


2008-04-24" 


580.9917 


-429.8 ± 


1.0 


451.1 


± 1.0 


2008-05-05" 


591.9748 


-428.9 ± 


1.0 


448.7 


± 1.0 


2008-05-07" 


593.9732 


-429.2 ± 


1.0 


448.0 


± 1.0 


2008-05-18" 


604.0442 


-429.6 ± 


1.0 


445.4 


± 1.0 


2009-01-07" 


838.1347 


-424.6 ± 


1.0 


381.3 


±0.9 



Notes, measurements from NACO have symbol ", whereas measure- 
ments from VISIR have v . The angles are all expressed in milliarcsec- 
onds. 

a different definition for the projection angle, leading to an ap- 
parent inconsistency with the other measurements. Transforming 
the Tango et al. projection angle PA using PA — > (180 - PA), 
these two data points become much more consistent with the 
other epochs and observing techniques. 

3.2. Orbital elements 

We adjusted the orbital parameters of the 5 Vel A-B pair to the 
whole sample of astrometric data, and the result is presented 
graphically in Fig. [8] The corresponding orbital elements are 
listed in Table [6] It should be noted that thanks to a semi-major 
axis twice more precise, and a period ten time more precise, the 
total mass value derived from based Kepler's third law is signifi- 
cantly improved, which becomes limited by our parallax estima- 
tion of jt = 39.8 + 0.4. 

M(Aab + B) = 6.15 ± 0.15 orb i, ± 0.17 parallax M (2) 



3.3. Physical properties of 6 Vel B 

We used a spectral energy distribution (hereafter SED) fit to 
the photometric data (Table [7} corrected for reddening assum- 
ing E(B — V) = 0.055. We used a carefully interpolated grid of 
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Fig. 8. left panel: Astrometric measurements of the relative position of r5Vel B (crosses) with respect to 5 Vel A (★ symbol), with 
our best fit orbital solution (solid ellipse, see Table [6] for the orbital parameters). Thin lines connect each measurement to the 
corresponding point on the adjusted orbit. The dotted line corresponds to periastron passage, and the diamonds corresponds to the 
position of 5 Vel B on Jan. 1 st every ten years between 1890 and 2020. Right pannels (labeled 1 to 4): corresponds to zooming boxes 
on the larger view. NACO and VISIR data (panel 1) are tabulated in Table |5J data in panel 2 comes from Horch, Franz & Ninkov 
(2000); data in panel 4 are from Tango et al. (119791 1. Zooming boxes' positions are reported on the main (left) panel. 



Table 6. Orbital parameters of the visual pair A-B. 



parameter 


this work 


Argyle et al. (2002) 


a CO 


1.996 ±0.012 


1.990 ±0.020 


e 


0.475 ± 0.003 


0.470 ± 0.020 


Period (yr) 


143.2 ± 1.2 


142 ± 13 


MJD 


51774 ±430 


51836.80 ±584 


i (deg) 


105.1 ±0.2 


105.20 ± 2.20 


£2 (deg) 


286.6 ± 0.36 


287.00 ± 1.30 


0) (deg) 


187.4 ± 0.6 


188.00 ± 14.00 


M(Aab + B) (M ) 


6.15 ±0.23 


5.88 ± 1.17 



Notes. The total mass (last line) is computed using Kepler's third law, 
assuming the revised parallax n = 39.8 ± 0.4 mas. 



Table 7. Spectral Energy Density (SED) fit of 6 Vel B, using our 
photometric measurements in the K and N bands (in two narrow 
band filters of VISIR, PAH1 and PAH2). 



filter 


meas. 


redd. 


SED 




modeled SED 








W/m 2 //jm 


W/m 2 /A«n 


V (Mag) 


5.54 


0.18 


2.159x10" 


10 


2.143 10- |U 


K(Mag) 


4.40 


0.02 


6.951x10" 


12 


7.135 10 12 


PAH1 (flux) 


0.94Jy 


negl. 


3.816x10" 


14 


3.816 10- 14 


PAH2 (flux) 


0.58Jy 


negl. 


1.369x10" 


14 


1.310 10~ 14 



Notes. The third column is the reddening correction corresponding to 
E(B - V) = 0.055. The model (last column) is for an angular diameter 
of 0.53 mas and effective temperature of 6600K. 



ATLAS models (e.g. Kurucz 120051) in order to derive the an- 
gular diameter and effective temperature of r5Vel B. We find 
a photospheric limb darkened angular diameter of 9ld(B) - 
0.530 + 0.011 mas and an effective temperature of T e ff(B) = 
6600 + 100 K. Based on our distance estimate, we can derive 
the physical radius to be R(B) = 1.43 ± 0.03 R Q and thus a lumi- 
nosity of L(B) = 3.5 + 0.2 L Q . Assuming the star is on the main 
sequence we can infer, based on the mass-luminosity relation by 
Torres, Andersen & Gimenez ( 120101) . that 6 Vel B has a spectral 
type F7.5V and a mass of: 

M photometr i c (B) = 1.35 ± 0.1 M B (3) 



These parameters estimated using an independent method are 
comparable to the values we obtained in Paper I. 

We can compare this mass estimate with the value we can 
compute from the mass of A (4.69 ± 0.03M o from table [2]) and 
A+B (6.15 ± O.23M from table[6j This leads to: 

M dynamic (B) = M(A + B) - M(A) = 1 .46 + 0.23 M Q (4) 

which is consistent with the photometric estimate (Eq.0 

4. Conclusion 

We presented a self consistent model of the triple stellar system 
6 Vel. Our model reproduces photometric, spectroscopic and in- 
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Fig. 7. c5 Vel Aa-Ab expected broadening functions computed for 
our models. For each star (Aa and Ab), the open dots are for the 
model (Table [3} and the thick line is for the analytical function 
resulting from the fit to the spectroscopic data (Fig. |2j. The small 
gray dots are for the model as well, but ignoring the gravity dark- 
ening. In the smallest panels are plotted the residuals "model" - 
"fit on spectroscopic data", with the dashed line being the best 
fit on the model using the same analytical function, showing the 
slight underestimation of the gravity darkening by our model. 
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terferometric data of the eclipsing pair Aa-Ab. We determined 
the orbital (Tabled and fundamental stellar parameters (TableO 
of the three components of the system. The physical properties 
of the eclipsing components are surprisingly similar to the A0V 
benchmark star Vega. Thanks to the resolution of the system us- 
ing the AMBER instrument, we also independently determine 
the distance to the system (n = 39.8 +0.4 mas, or 25.1 +0.25pc), 
as well as the interstellar reddening value towards this nearby 
system, with E(B - V) = 0.055. The combination of two fast 
rotating A stars of slightly different masses and a late F star, all 
coeval and with accurately measured fundamental parameters, 

will likely make of 5 Vel a cornerstone for the study of early- 

type main sequence stars. 2 Available at http:// www.jmmc.fr/amberdrs1 
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